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Abstract 


A three-dimensional fully viscous 
computer analyses, which retains the viscous 
nature of the Navier-Stokes equations, was 
evaluated to determine its usefullness in 
the design of supersonic inlets. This pro- 
cedure takes advantage of physical approxi- 
mations to limit the high computer time and 
storage associated with complete Navier- 
|q Stokes solutions. Computed results are pre- 
?! sented for a Mach 3.0 supersonic inlet with 
« bleed and a Mach 7,4 Hypersonic inlet. Good 
agreement was obtained between theory and 
data for both inlets. Results of a mesh 
sensitivity study are also shown. 
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total enthalpy 
metric coefficients 
mixing length 

reference length (35., SG cm 
(14 in.) for M3 inlet, 18.33 cm 
(7.22 in.) for P8 inlet) m 
bleed mass flow 
capture mass flow 
free-stream Mach number 
static pressure 
reference pressure (28.5 N/m 2 
(58.8 psf) for M3 inlet, 

701 N/m 2 (14.6 psf) for P8 inlet) 
Reynolds number per unit length 
based on free-stream (reference) 
conditions 

reference velocity (640.2 m/sec 
(2100 ft/sec) for M3 inlet, 

1222 m/sec (4008 ft/sec.) 
for P8 inlet) 

velocities in computational 
coordinates 

computational coordinates 
cartesian coordinates 
distance from wall along com- 
putational coordinate 
ratio of specific heats 
boundary layer thickness 
Von Karman constant, .43 
effective viscosity, laminar + 
turbulent 

turbulent viscosity 
density 
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Introduction 


The design of three-dimensional super- 
sonic inlets is a difficult task in view 
of the Wide operating range over which good 
performance is desired. Design of such an 
inlet system is strongly effected by the 
requirement that the aircraft operate nt 
all speeds from zero to the supersonic design 
point. In addition, most supersonic aircraft 
must have a subsonic cruise capability where 
fuel economy becomes important. As a result 
of these varied and sometimes conflicting 
requirements, the design of supersonic inlets 
is a difficult compromise. 

The foregoing considerations suggest 
that inlet design technology would benefit 
from a detailed and accurate flow field cal- 
culation procedure that includes shock- 
boundary layer Interaction affects. Within 
rectangular inlets operating at supersonic 
„ speeds, these viscous effects can bo classi- 
fied into three types: (1) incident shock 
reflections on the centerbody and cowl, 

(2) glancing shock interactions which take 
place along the inlet sidewalls, and (3) 
corner flow Interactions. Calculation pro- 
cedures which are either two-dimensional or 
do not account for boundary layer effects 
have limited application towards understanding 
these viscous interactions. 

A number of experimental and analytical 
papers have been published dealing with these 
shock interactions. The analysis and compari- 
son with data of an Incident shock wave 
reflection within an axisymmetric inlet have 
been accomplished by Fukuda, Hingst and 
Reshotkot and llingst and Towne 2 . In these 
analyses, control volume methods were used 
to solve for the properties downstream of 
the interaction zone given the conditions 
upstream of the interaction. A similar 
approach was used by Paynter’ for the analysis 
of the weak glancing sidewall interaction. 
These analyses yield correct trends and 
represent useful design techniques. An 
experimental study of the glancing sidewall 
interaction was done by Oskam, Vas and 
iiogdonoff4 and provides a detailed flow 
.’ascription from which to verify analyses. 
Except for the full Navier-Stokes analysis 
of the shock-corner boundary layer interaction 
by Hung and McCormack^ most of the analyses 
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of tills phenomenon ape Invlecid ! ' and provide 
no Information as to the viscous Interactions. 

There are however, a number of experimental 
reports 8 * 9 providing data on the viscous 
corner Interaction which are helpful 

understanding the flow structure itself. x - Momentum 


Although full Navicr-StokOB procedures 
would provide the necessary generality to 
predict the flow within three-dimensional 
supersonic inlets, the required computer time 
and storage would be prohibitive In terms of 
present computer technology. Combined viseoue- 
invisoid interaction analyses of the type 
suggested by Royhner and HiokcoxlO for axi- 
symmetrlc inlets could be implemented for 
three-dimensional inlets; however, the 
coupling procedures would bo very complex, 

An attractive alternative analysis would 
retain the general three-dimensional vie jus 
nature of the Navier-Stokes equations bu\ 
would take advantage of physical approxima- 
tions to limit the high computer run time 
and storage associated with a complete 
Navier-Stokes solution. In the present 
analysis, the assumption if made that there 
is a primary flow direction and that diffusion 
in this direction can he neglected. In this 
manner a set of equations are produced for 
fully viscous, predominantly supersonic, 
three-dimensional flow which can be solved 
by an efficient forward marching technique. 

This paper represents the first in a 
series of studies to evaluate the marching 
analysis, of Huggeln, Kroskovsky, and 
McDonald , designated PEPSIS, and to dete>r- 
mine its usefulness in the design of super- 
sonic inlets. Two inlets were chosen for 

this study; a Mach 3„0 configuration-*-* which 

obtains its compression by means of a crossed 
shock structure, and a Mach 7 .h hypersonic 
inlet- 1 --' 1 . Efforts to date hove concentrated 
primarily on the ramp/cowl shock wave boundary 
layer Interactions. 
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Governing Equations 

In this study, the inlet flow field is 
computed by a spatial marching method-*"*-*-*-* 1 ' 
which solves a simplified form of the three- 
dimensional Navier-Stokes equations. A 
curvilinear orthogonal coordinate system is 
used with coordinate directions x, y, and z 
and corresponding metric coefficients hj_, 

I12, and h 3 . Here x is defined as the 
streamwise or marching direction, and y and z 
are the two cross-flow directions. The equa- 
tions are first time-averaged so that they 
apply to both laminar and turbulent flow. 
Viscous terms are simplified by using an 
order-of-magnitude analysis. In particular, 
in order to allow the use of a marching pro- 
cedure, x-direction diffusion terms are neg- 
lected. The resulting equations are: 
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The energy equation is eliminated by 
assuming the total entholpy h 0 is constant. 
The pressure and density are related through 
the equation of state, 

p. ■~-f[h 0 --jp(u« + v*+w , l) 


A mixing length turbulence model is used, with 


The mixing length distribution is given by; 
i m (7) ■ 0.098 „ tort h[ <7 /(0.098 b )] 


The above aquations are solved by starting 
at the ramp tip and marching downstream 
using an alternatign-direction implicit 
technique. In this study, to avoid resolving 
the laminar sublayer, wall function boundary 
conditions were used to compute tangential 
volocities on all solid surfaces. Details 
of the solution procedure are published 
elsewhere 1 " 1 - 

In the above equations, the metric co- 
efficients can bn computed by any available 
method. For this study J since only inlets 
with rectangular cross-sections and parallel 
sidewalls are considered, h^ and hg are 
functions of x and y only, and 113 * 1. The 
mesh was therefore generated using the two- 
dimensional method of Anderson 1 -?. 

Test Cases 


In order to evaluate the computational 
procedure described, two test cases were run 
with the PEPSIS analysis. The first configu- 
ration considered was a large scale, variable 
geometry, mixed compression inlet * 2 designed 
for Mach 3.0 and a Reynolds number per meter 
of 7.2X10° (2.2X10° per foot). This test- 
configuration was designated the M3 inlet 
configuration. It was chosen to study the 
behavior of the PEPSIS analysis under condi- 
tions of high Reynolds number and with a 
complex crossed shock wave pattern in a 
largely inviscid core. The second case was 
a Mach 7.4 hypersonic inlet 1 - 3 designated 
the P 8 inlet. This case was chosen to study 
the behavior of the PEPSIS analysis under 
conditions of strong shock wave interactions 
with thick boundary layers. In addition, 
this inlet was used by Knight 1 -® as a test 
case for a two-dimensional Navier-Stokes 
solution. Each inlet configuration was run 
with two computational mesh distributions, 
Table I, designated the coarse and medium 
mesh solutions. The coarse mesh had 20 grid 


points distributed between the ramp and cowl 
surfaces while the medium mesh had 45 grid 
points. The mesh points for both the coarse 
and medium mesh solutions were packed In the 
region near the ramp and cowl surfaces to 
resolve the wall boundary layers. In this 
study, since only 6 mesh points were used 
between the inlet sidewall and the symmetry 
plane, the sidewall boundary layers were 
not resolved. This did not lead to any 
computational problems, however. For both 
inlets, the streamwise step size distribution 
was the same for the coarse and medium mesh, 
hater steps in the evaluation procedure will 
study the effects of varying the spanwlee 
and streamwise mesh size. The total computing 
times, i.o., CPU plus I/O time, on the UXIVAC 
1100/42 are also listed in Table I. 

M3 Inlet Configuration 

A schematic diagram of the M3 inlet con- 
figuration allowing the inviscid shock 
structure is presented in figure (1) . The 
inlet capture area was 1264 cm 2 (190 in 2 ) and 
the initial ramp was Inclined 7 degrees to 
the horizontal. The second ramp was inclined 
14 degrees to the horizontal while the cowl 
surface was parallel to the horizontal. The 
inlet achieves internal compression by 
means of a crossed shock structure as indi- 
cated in figure (1) . Calculations were made 
on the M3 inlet with the ramp, cowl and 
sidewall bleed regiong sljpwn in. figure ( 2 ). , 

The bleed rates were computed from the data 
for configuration 80 of reference 12 and are 
listed in Table II, A comparison between 
the coarse and medium meshes used in the 
calculation of the M3 inlet configuration 
is shown in figure (3) . A section upstream 
on the romp tip was added to properly resolve 
the Initial ramp shook in the calculations. 

In the description that follows, u and v are 
the streamwise and cross-flow velocities 
in this coordinate system. In addition, the 
gapwise distance is defined as the distance 
along a coordinate line from the ramp to 
the cowl. Comparisons between the coarse 
and medium mesh solutions are presented in 
figures (4)-(0). In figure (4), computed 
Mach number profiles in the center plane of 
the M3 inlet have been plotted within the 
inlet geometry. The coarse mesh solution 
exhibited substantial numerical noise, par- 
ticularly in those regions where an oblique 
shock either formed or was reflected from a 
solid surface. This can be seen more clearly 
in figure (5) where the u-velocity (stream- 
wise velocity) profiles have been plotted 
as a function of normalized gap distance 
across the flow field. Differences in the 
computed shock wave structure are also evident- 
in figure (5) . A more distinct second ramp 
shock is apparent in the medium mesh solution 
in addition to an upstream shift in the first 
ramp shock-boundary layer interaction. The 
dramatic differences between the coarse and 


medium mesh solutions con also be seen in 
figure (0) which presents the v-velocity 
(normal velocity) profiles in the M3 inlet 
center plane as a function of normalized 
gap distance. The v-velocity is very sensi- 
tive to mesh resolution and it is evident 
from figure (6) that the numerical problems 
encountered in the coarse mesh solution were 
essentially eliminated in the medium mesh 
solutions. 

A comparison between the two solutions 
and experimental static pressure measurements 
is presented in figure (V) . The medium mesh 
gave a forward shift in the shock-boundary 
layer interactions, but both solutions agree 
well with the data. The boundary layers on 
both the ramp and cowl surface were thin 
relative to the channel gap dimensions so 
that detailed resolution of these interactions 
would require additional mesh points, The 
solution ended just upstream of the throat 
because of a substantia', region of subsonic 
flow in the core of the flow field which the 
PEPSIS analysis cannot presently handle. 

In order to continue the solution past 
this region, a second medium mesh calculation 
was performed at over speed conditions, i.e., 
at a free stream Mach number of 3.25, and 
the results are presented in figure (8) . 

The Mach number profiles which are plotted 
within the inlet geometry reveal that the shock 
impinging on the romp shoulder has been can- 
celled by the sudden expansion around the 
shoulder. Compression downstream of the 
shoulder thus takes place by a single reflect- 
ed shock system. The increase in boundary 
layer thickness as a result of the sudden 
expansion around the ramp shoulder is also 
evident. Again to properly resolve the flow 
phenomena in the region of the ramp shoulder 
would require a much greater ir.jsh density. 
However, the PEPSIS analysis does simulate 
the gross features of the flow with the pre- 
sent mesh resolution which makes it ideal 
for design. 

P8 Inlet Configuration 

In order to investig. ste the ability of 
the PEPSIS analysis to compute the details 
of tbs shook wave-boundary layer interactions 
a series of calculations was made on a Mach 
7. 4 hypersonic inlet. Under hypersonic condi- 
tions, much thicker boundary layers are 
encountered so that greater resolution of 
these interactions can be achieved without 
increasing the number of mesh points. 

A schematic diagram of the P8 hypersonic 
inlet showing the inviscid shook structure 
is presented in figure (9) . The initial ramp 
angle was inclined 6.45 degrees to the 
horizontal while the initial cowl angle was 
inclined -1 degree to the horizontal. The 
leading edge of the cowl on the P8 inlet was 


blunt, with a diameter of 0.114 cm (.0449 in). 
Since it is impractical to resolve this 
blunt body flow in a complete inlet analysis 
an effective sharp cowl leading edge was 
defined. Using the expressions from Hayes 
and Probstciiv LJ for two-dimensional shock 
standoff distance, and positioning the 
initial cowl shock tangent to the bow lip 
shock, figure (10), gave an effective dis- 
placement of 0.26 cm (.10 in). This compares 
to an effective cowl lip displacement of 2.B 
cm (1.1 in) used in the two-dimensional 
Navier-Stokes solution of this inlet by 
Knight^. 

A comparison between the coarse and 
medium mesh Mach number profiles plotted 
within the inlet geometry is presented in 
figures (11) and (12) . Figure (13) shows 
the developing Mach number profiles for 
the entire P8 inlet, while figure (12) pre- 
sents an enlarged view internal to the inlet. 
More distinct initial ramp and cowl shock 
waves are apparent in the medium mesh solu- 
tion. With the enlarged view of the internal 
Mach number profile development, figure (12) , 
the tjwl shock wave interaction with the 
thick ramp boundary layer becomes evident. 

The reflected shocks off both the ramp and 
cowl boundary layers are also clearly pre- 
sent. Differences between the coarse and 
medium mesh solutions become more apparent 
by plotting the Mach number profile develop- 
ment in terms of normalized gap distance 
across the inlet channel, figure (13) . , 

The numerical noise disturbances generated 
just downstream of both the ramp and cowl 
shock waves in the coarse mesh solution, 
figure (13a) , are propagated along Mach 
lines and act like weak shock waves. 

Increasing the mesh density essentially 
eliminated these numerical disturbances. 

The medium mesh solution, figure (13b) , 
had enough grid points to resolve the detailed 
flow events within the length scale associated 
with this interaction. The boundary layer 
buildup along the ramp surface and its 
interaction with the cowl shock is clearly 
evident. The static pressure profile deve- 
lopment in the center plane of the P8 inlet, 
figure (14) , reveals more clearly that the 
numerical noise disturbances are propagated 
along Mach lines in the coarse mesh solution 
and that increased mesh resolution eliminates 
these problems. Note also that there is a 
substantial normal pressure gradient across 
the ramp and cowl boundary layers within the 
interaction zone. 

Comparisons of the coarse and medium 
mesh static pressure solutions with measure- 
ments in the center plane of the P8 inlet 
are shown in figure (15) . The calculations 
were performed assuming turbulent boundary 
layer flow along both the ramp and cowl 
surfaces while the experimental test had 
substantial laminar flow on both surfaces. 


As a result, the analysis predicted a thicker 
boundary layer entering both Interaction 
regions which resulted In a slight upstream 
shift In the shock location. In spite of 
these differences very good static pressure 
predictions were achl- ved with the PEPSIS 
analysis. The PEPSIS code is presently being 
updated to Include a laminar-turbulent 
transition capability, it is expected that 
with thinner boundary layers still better 
agreement will ho achieved. A comparison 
between the PEPSIS viscous analysis and 
invlseld calculations 2 ') Is presented In 
figure (10) and shows a strong ramp shock 
Interaction extending well upstream of the 
invlseld shock impingement point and a 
substantial upstream displacement In the 
cowl shock interaction region. 

Concluding Remarks 

A scries of calculations were performed 
on two inlet configurations using a three- 
dimensional viscous flow analysis (PEPSI8) 
and the following remarks can be made: 

1, With proper overall mesh resolution, 
the PCPSIS analysis is able to define the 
basic features of high Reynolds number thin 
boundary layer flows with highly complex 
shock wave structures in a largely invlseld 
core. 

2, Even when the solution mesh is too 
coarse to properly resolve the small length 
scale events associated with thin boundary 
layer-shock wave interactions, the overall 
features of these evento can be predicted. 

3, With adequate mesh resolution of the 
shock wave-boundary layer interactions, the 
PEPSIS analysis predicts the gross physics 
of this flow event, 

4, Depending on the mesh density, the 
PEPSIS analysis can be used either as a 
design tool, where the overall inlet flow 
field is the important factor, or as an 
analysis tool where the small scale flow 
events ocouring in inlets can be studied. 

5, Numerical noise disturbances can 
develop within the solution which will 
distort the flow field, but proper mesh 
resolution can eliminate this problem. 
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TABLE I. - TEST CASES FOR PEPSIS 


Test cases 

Computational mesh 

Computing 
time* 
on UNIVAC 
1100/42, 
min 


Distribu- 

tion 

Total 

M3 inlet 




Coarse 

20X6X180 

21 600 

58 

Medium 

45X6X177 

47 790 

129 

Medium 

45X6X260 

70 200 

190 

P8 inlet 




Coarse 

20X6X350 

42 000 

113 

Medium 

45X6X340 

91 800 

248 


*CPU plus I/O time, 


t * 

TABLE II. - BLEED RATES FOR M3 


INLET CONFIGURATION 


Bleed 

zone 

Bleed mass flow, 

V*- 

Bleed velocity, 
v/u re£ 

r a 

0.0097 

0.00902 

r b 

.0053 

.00528 

R C 

.0193 

.00769 

R D 

.0214 

.00832 

C A 

.0120 

.00295 

C B 

.0699 

.00431 

S A 

.0120 

.00451 

S B 

.0147 

.00748 

s c 

.0176 

.01277 


m = 7.006 kg/sec (15.45 lb /sec). 

oo m 

u re £ = 640.2 m/sec (2100 ft/sec). 
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Figure 4 - Mach number profiles In center plane of M3 .Inlet configuration. Mo, ■ 3,0; Re • 
7.2xl0 6 /m <2. 2xl0 6 /ft). 
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ment was obtained between theory and data for both inlets. Results of a mesh sensitivity study 
are also shown. 


17. Key Words (Suggested by Author(s)) 

Supersonic inlet 
Parabolized Navier-Stokes 
Shock waves 

Viscous flow 
Marching analysis 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 34 

19. Security Classif, (of this report) 
Unclassified 

20, Security Classif. (of this page) 

Unclassified 

21, No. of Pages 

22. Price' 


' For sale by the National Technical Information Service, Springfield, Virginia 22161 


